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Abstract 

We consider the time-dependent electron transport through a quantum dot cou- 
pled to two leads in the presence of the additional over-dot (bridge) tunneling chan- 
nel. By using the evolution operator method together with the wide-band limit 
approximation we derived the analytical formulaes for the quantum dot charge and 
current flowing in the system. The influence of the external microwave field on the 
time-average quantum dot charge, the current and the derivatives of the average 
current with respect to the gate and source-drain voltages has been investigated for 
a wide range of parameters. 

1 Introduction 

Electronic transport in mesoscopic systems has been at the focus of experimental and 

theoretical interest during the last decade due to recent development in fabrication of 

small electronic devices and their interesting equilibrium and non-equilibrium properties. 

Especially interesting are the transport properties of a quantum dot (QD) under the 

influence of external time-dependent fields. The high-frequency signals may be applied to 

a QD and the time-dependent fields will modify the tunneling current. 

New effects have been observed and theoretically described, e.g. photon-assisted 

tunneling through small quantum dots with well-resolved discrete energy states [1, 2, 

3], photon-electron pumps [4, 5, 6] and others. One can investigate the current flowing 
"'^corresponding author, e-mail: taranko@tytan.umcs.lublin.pl 



through a QD under periodic modulation of the QD electronic structure [7] or periodic 
(non-periodic) modulation of the tunneling barriers [6] and electron energy levels in both 
(left and right) electron reservoirs [8] (see also [9, 10]). 

The progress of nanomaterials science has enabled the experimental study of the 
phase coherence of the charge carriers in many mesoscopic systems. The asymmetric 
Fano line shapes [11] are observed whenever resonant and nonresonant scattering paths 
interfere. In some nanostructures, e.g. in single-electron transistors, the Fano resonances 
in the conductance were observed [12], which imply that there are two paths for transfer of 
electrons between a source and a drain. Especially the recent experimental and theoretical 
study with a low-temperature scanning tunneling microscope (STM) of the single magnetic 
atom deposited on a metallic surface showed the asymmetric Fano resonances in the 
tunneling spectra [12-15]. The STM measurements indicate that in tunneling of electrons 
between STM tip and a surface with a single impurity atom two different paths are present. 
The electrons can tunnel between the tip and the adsorbate state and directly between 
the tip and the metal surface. The electronic transport through a QD coupled to the 
electron reservoirs within a model with two electron tunneling channels was considered in 
Ref. [16] and it was shown, that transport of electrons through both channels leads to an 
asymmetric shape of the zero bias voltage conductance curves, which is typical behaviour 
for a Fano resonance resulting from constructive and destructive interference processes 
for electrons transmitted through both channels. 

In all papers mentioned above and relating to the electron transport through a QD 
with the additional (bridge) transmission channel the external fields were not applied 
and the considered systems were driven out of equilibrium only by means of a dc volt- 
age bias. In this paper we address the issue of a QD with a bridge channel between a 
source and a drain driven out of equilibrium by means of a dc voltage bias and additional 
time-dependent external fields. In this manner, our paper can be seen as generaliza- 
tion of Ref. [9] to the case of a QD with the additional bridge channel in the presence of 
external microwave fields which are applied to the dot and two leads, respectively. In liter- 
ature, different theoretical approaches have been developed to treat the time-dependent, 
nonequilibrium electron transport processes in the mesoscopic systems. It seems, that 
the most popular is the nonequilibrium Green's function method. However, these Green's 
functions depend on the time arguments and for non trivial quantum models it is a rather 
difficult task to calculate them. In our treatment of the time-dependent tunneling through 



a mesoscopic system we use the evolution operator technique (e.g. [17, 18]). The final 
expressions for the QD charge and the current flowing in the system can be described 
in terms of the corresponding matrix elements of the evolution operator. In our earlier 
work [19] we have considered the similar problem solving numerically the corresponding 
sets of the differential equations satisfied by the matrix elements of the evolution oper- 
ator. Due to the complexity and a large number of these equations we considered only 
a very limited number of intereting cases although we were able to take into considera- 
tion the electronic structure of the lead energy bands and the specific time-dependence 
of the QD-lead barriers. Here we give the analytical expressions for the QD charge and 
current assuming so called the wide-band limit approximation and the time-independent 
strength of the QD-lead barriers. As a result, due to these analytical forms, we are able 
to analyze the required characteristics of the considered system for the very broad class 
of parameters. Additionally, due to the final forms given for some matrix elements of the 
evolution operator it is possible to build up the expressions for the QD charge or current 
in the form of the perturbation series. 

In the next Section we present the model and formalism and give the resulting 
expressions for equations for the corresponding matrix elements of the evolution operator. 
In Section 3 we obtain the approximate solutions for all required matrix elements and give 
the final forms for the QD charge and the current flowing in the system. The results are 
presented in Section 4 which includes also the summary and a brief discussion. 

2 Model and calculation method 

We model the QD coupled to the left and right electron reservoirs with the additional 
bridge tunneling channel between them by the usually used Hamiltonian H = Hi + V, 
where 

^1 = Jl^kJ^)(^t^kc. + £dit)a^ad , (1) 

ka 

^ = E VUtK<^^ + h-c- + E VhknitK^k, + h-c. (2) 

The operators {a^ ), ad{a^) are the annihilation (creation) operators of the electron 
in the lead a {a — R) and in the QD, respectively. The couplings between QD and 
lead states and between both lead states are denoted by Vr , and Vr t , respectively. For 



simplicity, the dot is characterized only by a single level Ed and we have neglected the 
intradot electron-electron Coulomb interaction. We assume the case in which there exist 
microwave fields applied to the leads and QD. In the adiabatic approximation our time- 
dependent driven system is described by Sji^if) = e^^ + Aq, cos cut, Edit) = Ed + AdCosa;t, 
i.e. the energy levels of the leads and QD are driven by the ac field with the frequency uj 
and the amplitudes and A(i, respectively. 

We describe the dynamical evolution of the charge locahsed on the QD and the 
current fiowing in the system in terms of the time evolution operator U {t, to) (in the 
interaction representation) which satisfies the equation 

t^U{t,to)^V{t)U{t,to), (3) 

where 

Vit)^Uoit,to)Vit)U+it,to), (4) 

Uo{t)^Texp (^lJ dt'Hi{t')j . (5) 

Here we assume that the interaction between QD and leads and between both leads is 
switched on in the distant past to, i.e. V^^^{t) and ^^jfc^(^) equal to zero for t < to and 
takes constant values for t > to- 

The QD charge and currents fiowing in the system can be obtained from the knowl- 
edge of the appropriate matrix elements of the evolution operator U {t, to)- The QD charge 
is given as follows (cf. [18]): 

ridit) = nd{to)\Udd{t,to)\' + i:nj:Jto)\U,j;Jt,to)\' , (6) 

where Udd{t-,to) = {d\U{t,to)\d) and U^^^{t,to) = {d\U{t,to)\ka) denote the matrix el- 
ements of U {t, to) calculated within the basis functions containing the electron single- 
particle functions of the leads and QD, Iki), and \d), respectively. nd{to) and n^^ito) 
represent the initial filling of the corresponding single-particle states. 

The tunneling current fiowing, e.g. from the left lead into the QD and the right 
lead, ihit) , can be obtained from the time derivative of the total number of electrons in 
the left lead, j^it) = —ednL{t)/dt (cf. [9]), where 

Mt) = E^fc,W = EM^o)|C/fc,,(i,io)r 

+ E%.(^o)it/^„,-.,(^,io)r+E^.-«(^o)it/^„,-jt,to)r). (7) 

9i kR 



Let us begin with the calculations of the QD charge nait). Then we have to calculate 
the matrix elements Udd{t,to) and U^j^^{t,to). Using the identity operator I — \d){d\ + 

— * — * 

\ka){ka\ the following set of coupled equations can be obtained from Eq. (3): 

^Ud^it, to) = E ViJt)U^,Jt, to) , (8) 

§iUkJt^ ^o) = -iV^Jt) U^it, io) - ^ E \~,,m-,Jt, to) , (9) 
^U^kjt^to) = -iV^Jt) Udd{t,to) -ij:\R,{t)Uj;Jt,to) , (10) 

where 

Vkjt) ^ {K\V{t)\d) = V^,^,exp{t{e^-ed){t-to) 

—i{ Ad — Aa) {sin cut — sin u!to)/uj) , (11) 

^ {kL\Vit)\kR)^V^j^exme^,^-ej^Jit-to)- 

— {{Al — Aji){sin Lut — sin u!to)/u!) . (12) 

It is easy to show that the equation for Udd{t, to) can be written as follows: 
dUdd{t,to) 



^-j dt'mt, t')Udd(t', to) + Y:jCi;Jt, t') Uj;Jt', to)) , (13) 
where 

^M-i:vijt)v,jt'), (14) 

^a;.(^>0 = Ev^,(^)%.^,(0, (15) 

and the similar equation can be written for ^^^(i, ^')- 

The formal solution of Eq. (9) written in the form 

t t 

UUt^ to) = -tJ dt%Jt')Udd{t', ^o) - ^ E / dt'\j;^{t')Uj;Jt') , (16) 
to kR to 

(and the similar equation for U^^^{t,to)) can be iterated giving 
t 

ULd{t,to) = -i j dtiVL{ti)Udd{tuto) 

to 

+ E(-^)' E / "^^^ • • • / ^^1^ ^^1^2 ^^2) ■ ■ ■ K,.^ {tj) ■ Udditj, to) . 

j=2 i?i,L2,i?3,i4,---,aj-ltQ to 

(17) 



Here we have introduced the abbreviated form for the vector and replaced it by a 
and Vi^^^{t) = Va{t)- The equation for Ujid{t,to) can be obtained from Eq. (17) by 
interchanging L R. Inserting these expressions for ULd{t-ito) and URditjto) into Eq. 
(13) one obtains the closed integro-differential equation for Udd{t,to) 

d 



Q-^Udd{t,to) = -J dhlC{t,h)Udd{tuto) 



to 

oo * *Jt1 

- E Jdh-.J dt,Vl^{t)VL,R,{h)...V^,{t,)Ud4{t,,to) 

3=2 Li,R2,L3,--;ajtQ to 

+ (the second term with the change L ^ R). (18) 

Equation (18) together with the expressions for Vl(^) and VLR{t), Eqs. (11, 12), and for 

)C{t,ti) written as follows 

lC{tM)^ E |Kr2^a(i-^i) exp(i£d(i-ii)+i(Ad- A,,)(sina;i-sina;ti)/u;) (19) 

a=L,R 

where Va{t) a Fourier transform of the cc-th lead density of states, gives exact, closed 
equation for Uddit,to)- Here we have assumed that Vj^^ does not depend on the wave 
vector ka and then the similar assumption will be made for V^^^^. 

Under the wide-band limit (WBL) approximation (e.g. [9]) this equation can be 
analytically solved and such solutions will be considered later. Formally, solving Eq. (18) 
and inserting its solution to Eq. (17), the solutions for ULd{t,to) and Ujid{t,to) can be 
obtained. These functions are needed in the calculations of the first term oijL{t) (see Eq. 
7). 

In order to calculate nd{t) we still need Uda{t,to)- Writing down the set of closed 
equations for UdR{t,to), Uji^ji^{t,to) and ULR{t,to) (obtained on the basis of Eq. (3)) and 
performing similar calculations to those described above, one obtains for UdR{t, to) (and 
similar equation for UdL{t,to) by interchanging L ^ R): 

d , . 

-Q^UdR[t, to) — 

t 

- iV*(t) + (-i)2 / dt,^Vl{t)VLR{t,) 

to ^ 

OO J ^ 

+ E i-^y^' E / • • • / dt,Vl{t)Vn,L,{t,) . . . Vr.^n{t,) 

j=2,4,6,... Ri,L2,R3,-,Lj^^ 

CO 1^ J ^ 

+ E Hy^' E Jdti--- 1 dtj Kit)VL,RAti) ■ ■ ■ VlM^j) 

j=3,5,7,... Li,R2,L3,...,Lj 



i ^ 

+ Ei-^y^' E Jdti--- I dtj Kit)^LMti) ■ ■ ■ VaMUaR{t,, to) , (20) 

j=2 Li,R2,L3,...,Q!jJq 

and for ULR{t,to) (needed for calculation of the last term of niit), Eq. (7)): 
t 

ULR{t,to) =-l j dhVLRih) 



to 

ti-i 



+ E Hy E Jdh jdh... J dt^vUhWR.M ■ ■ ■ VL^-Mtj) 

j=3,5,7... Ri,L2,R3,--;Lj-n^ 
t 

- i I dtiVL{ti)UdR{ti,to) (21) 

to 

+ E(-^)' E J dtij dt2... J dtjVLR,{ti)VR,L,{t2)...K,{tjW^ 

j=2 _Ri,L2,.--,Q!j-lio to to 

The analytical solutions of these equations under the WBL approximation will be dis- 
cussed in the next section. 

For calculation of n-^^ (t) one needs the functions U Ld{t, to) , U LR{t, to) and UL^L^{t,to). 
The first two functions are given in Eqs. (17) and (21) and UL^L^{t,to) should be calcu- 
lated from the set of coupled equations for UdL{t,to), ^/liL2(^)^o) and URL{t,to) obtained 
from Eq. (3). The result is as follows: 

oo * *l 

'it,to) = Sl,L + E i-^y E jdtijdt2... I dtjVL,RAtl)VR,L,{t2) . . . VR._^L(tj) 

J=2,4,... Rl,L2,R3,:.,R3--Lto to to 



t 

-I 

to 



i I dt,Vl^{t,)UdL{ti,to) (22) 



+ E(-^)' E J dtij dt2... J dtjVLMtl)yR.L,{t2) . . . Va,_,L{tj)UdL{tj, to) , 

3=2 Ri,L2,R3,-,aj-n^ 

where UdL{t,to) is given by solving Eq. (20) with the replacement L ^ R. 



3 Analytical solutions in the WBL approximation 

In order to calculate the QD charge nd{t) or current jhit) one has to solve, in the first 
step, the integro-differential equations satisfied by Udd{t,to), Eq. (18), and UdR/L{t,to), 
Eq. (20). The other needed functions ULd{t-,to), UL^Liit^to) and ULR{t,to) can be ob- 
tained from Eqs. (17, 21, 22) inserting into them Udd{t,to), UdL/R{t,to) and performing 



multiple time and A;- vector integrations. Unfortunately, it is a very difficult task to solve 
the integro-differential equations and to perform these integrals in a general case when the 
leads are characterized by some density of state curves. Here, we use the WBL approxi- 

— * 

mation, under which all multiple time and fc-vectors integrals can be performed without 
difficulties. The WBL approximation has been widely used in calculating of many prop- 
erties of mesoscopic systems (e.g. [2, 3, 5, 9, 10]). It is justified under the conditions that 
the QD energy level linewidth is much smaller than the bandwidth of the leads and that 
the density of states and hopping matrix elements vary slowly with energy. Furthermore, 
as we are not going to consider the case of the QD energy level lying close to the edges of 
the leads energy band or lying close to some singular structure present in the leads den- 
sity of states, then application of the WBL approximation should be fully justified. The 
conditions under which we perform our calculations are satisfied in most experimental 
constructions of mezoscopic systems. As a check, we have performed the direct but time 
consuming numerical integration of Eqs. (8-10) (and similar equations for other func- 
tions) for the rectangular leads density of states and did not find any differences in the 
results for the time-averaged QD charge or the currents fiowing in the considered system. 
Therefore, in this section we consider the electron transport through the QD with the 
additional bridge tunnel over the dot within the WBL approximation. In this approach 
the solutions of the intcgro-diffcrcntial Eqs. (19) and (20) in the analytical form can be 
obtained and the infinite sums of all terms in Eqs. (17) and (21) can be calculated. 

Let us consider the equation for the matrix element Udd{t, to), Eq. (18). The function 
K,{t,t'), Eq. (19), is approximated in WBL as follows 

/oo 
deVccie) exp(-i£(i - ti)) exjp{-ied{t - ti) 
-°° 

-\-i{/S.d — Aa)(sina;i — sva.ujti)/uj) =^ 



^ E^27r5(t-t0, (23) 

where the leads density of states 'Da{e) was replaced by the rectangular density of states 
with the "effective" bandwidth 

Using similar approximation in calculations of the multiple integrals present in Eq. 
(18) one obtains 



-Q^ddit, to) = \^-- - J^i-'^^RL/Dy I Udd{t, to) , (24) 

where r„ = 27rVyD, T = Tl + Tr, = Dr = D, V = V^. 



Assuming ttVrl/D < 1 the series can be summed up and finally the equation for 
Udd{t,to) reads 

^Uddit, to) = -CiUddit, to) , (25) 

where Ci = {27rV^/D)/{l + inVRL/D). 
Similarly, Eq. (20) becomes 

^Udait, to) = -iVdd{t)/{l + iVrl'k/D) - C,Uda{t, to) . (26) 

The solutions of Eqs. (25) and (26) are as follows: 

Udd{t,to)=exp{-Ci{t-to)), (27) 

t 

Uda(t, to) = + iVrltt/D) I dtiVdaih) exp(-Ci(i - ti)) (28) 

to 

and the QD charge nd{t) can be easily obtained, Eq. (6). It can be verified, that the first 
term of Eq. (6) tends to zero as t — to ^ oo and finally for the QD charge we have 

where 

t 

Aa{e,t) — — j dti ex.p(i{ed — s){t — ti) — i{Ad — A.a){smu;t — smu;ti)/u; — Ci{t — ti)) (30) 
to 

In the limit of vanishing bridge over the QD, Eq. (29) reproduces the result of Ref. 

[9]. 

The current Jl(^) flowing from the left lead into the QD and the right lead is calcu- 
lated from the evolution of the total number of electrons in the left lead (see Eq. 7) and 
one can read: 

jr.{t) = '^ReLdito)Y.UlJt,to)jU^,Jt,to) 
+ E n,-,(to)C/i,,-.,(i,to)|c/g„,-,(^,^o) 



+ E %^.(^o)t/|,,,^^(t,to)|t/,-,,^,(t,to) 



(31) 



The functions Uj^^^, ^k^qL ^kLk^i^^^o) are calculated according to Eqs. (17), (22) 
and (21), respectively, and after summing up of the corresponding multiple integrals one 



obtains: 

t 



UkJt,to) = -Y^JdtiVj:Jh)U,,{h,to), (32) 

to 

to 

t 

J dtiVkJh)Ud,-,(h,to), (33) 



1 + ix 



to 

t 



to 
t 

to 

where Udd{t,to), U^^^{t,to) and U^^Jt.to) are given in Eqs. (27, 28) and x = Wrl/D. 
One can verify, that the first term in Eq. (31) tends to zero ior t — ^ oo as we have 
for this term 

n,(to)|l/(l + ix)\^ r/2 exp(-r(t - to) ReC^) , (35) 

where ReC2 = Re{l/{1 + ix)) > 0. 

The second and third terms of the expression for jhit), Eq. (30), together with Eqs. 

(27, 28, 31-33) give finally the time averaged current: 

1 2t^ c r/2 

In the vanishing bridge channel case Eq. (35) coincides with the results of Ref. [9]: 

{jl'''='it)) = ^{n,{t)) + de h{e)lm{A,{e,t)) . (37) 

The current {jL{t)), Eq. (36), fiowing from the left lead to the central region and to the 
right lead (through the bridge channel) is the superposition of four terms. The first term 
corresponds to the current between two leads and this term is not disturbed by the QD. 
The form of the second and third terms is the same as for (jr^^=°(t)), Eq. (37), except for 
the renormalization constants due to the additional tunneling channel. Note, that some 
additional renormalizations also occur due to Vrl which enters into the expression for 
nd{t), Eq. (29), and for Aa{£,t), Eq. (30). The last term of Eq. (36) is the interference 
term due to the simultaneous tunneling through two channels. 



4 Results and discussion 



Here we show the numerical results of the time-averaged QD charge {nd{t)) and the current 
{jhif)) and its derivatives with respect to the QD energy level position and the chemical 
potential (or equivalently, with respect to the gate and source-drain voltages) for 
different sets of parameters which characterize our system. We assume the temperature 
T = K and the time-average of time-dependent quantities fit) is defined by 



and because f{t) is a periodic function of time, wc average it over the period 2ti/uj. We 
take the chemical potential of the right lead as the energy measurement point, ^r — Q. 
As the potential drop between the left and right leads is given hy /il — /J-r — eV^-d and 
Vs-d is the measured voltage between a source and a drain, then the derivatives of the 
current with respect to fi^ will correspond to the derivatives d{jL{t)) /dVs-d usually 

measured in experiments. In experiments the gate voltage controls the position of the QD 
energy level Sd (regardless how complicated the relation between the gate voltage and Sd 
is) and for that reason to mimic measurements of the QD charge or current vs. the gate 
voltage we have calculated them vs. the position of the QD energy level Ed- 

The values of the hybridization matrix elements ^ present in the Hamiltonian do 
not enter the final expressions for the current or QD charge obtained within the WBL 
approximation. Usually the effective linewidth = 27rE^^|Vg^^|^(5(£J — £g^) is introduced. 
However, in our calculations the others hybridization matrix elements appear, V^^ -^^ = 
Vrl, responsible for the additional tunneling channel for which we should take some values 
in order to perform numerical calculations. We have taken the values comparable with 
Vj^^^ and estimated Vg^^ (assuming its /c-independence, Vg^^ = Va = V) using the relation 
Fq, = 27r| Vq,P/-Dq, where Dq, is the a-lead's bandwidth and — 100 Fq, (F^ — Tr — T, 
Dl — Dr — D was assumed). We assumed the amplitude of the QD energy levels 
oscillation to be one half of and Ar — 0, if otherwise stated. In our calculations 
the values Vrl were taken from the range (0-10), in F units. 

In the first three figures. Figs. 1-3, we present the overall shape of the average 
current (jl(^)) and the derivatves of the average current with respect to the QD energy 
level Ed and the left chemical potential /il, respectively, against Ed and /il. The upper 
panels correspond to the Vrl — case and the lower ones present the results obtained for 



T 




(38) 



— T 



the non- vanishing over-dot tunnehng channel, Vrl 7^ 0. There are quite visible differences 
between the case of a QD with and without additional bridge tunneling channel (cf. upper 
and lower panels in Figs 1-3) . For better visuahzation of the peculiarities of the presented 
functions and for simpler discussion let us consider the specific cuts of the surfaces given 
in Figs. 1-3. 

At first, let us consider the dependence of vs. Sa for given values of the left 

lead chemical potential ii^. In Fig. 4 we present such curves for different values oi /il- the 
subsequent curves beginning from the lower one correspond to /il — —4 up to the upper 
curve (with the step Afii = 1) obtained for /xl — 8. The left (right) panel corresponds 
to Vrl = {Vrl = 10). In the case of vanishing over-dot tunneling channel (the left 
panel) the current has a simple structure - a single peak localized in the middle between 
/Il and /ir for smaller values oi /il. The width of this peak increases with increasing |//^| 
{iir — 0) and for greater values of the current is almost independent of £d localized 
inside the energy region between fiji and fii- For non-vanishing over-dot tunneling (Fig. 
4, the right panel) the curves (jl(^)) become asymmetric. With increasing source-drain 
bias, the current possesses greater values in comparison with the Vrl — case due to the 
direct tunnehng between both leads. Note, however , that due to the interference effects 
the resulting {jL{t)) curves are asymmetric. The interference effects are most visible 
approximately for Sd lying in the region {fiR,fiL). 

In Fig. 5 we show the average current {jhit)) vs. the left lead chemical potential 
for several values of Sa- For vanishing Vrl the corresponding curves are nearly asymmetric 
with the asymmetry point /il — s:d- With increasing //^ at fixed ea the current achieves 
a constant value depending on the position of the QD energy level ea with respect to the 
Hr = 0. It means, that electrons which occupy the lead energy levels not too distant from 
Ed take part in the tunneling process. For greater fi^ most of the lead energy levels lying 
far away from the Ed are inactive in the tunneling between leads through the QD energy 
level. However, at non-vanishing Vrl (see the right panel of Fig. 5) these lead energy 
levels are active and the current (jL(t)) vs. is of much richer structure. The current 
is nearly linearly growing with the increasing fiL (for larger fii) because the tunneling 
through the QD can be neglected compared with the direct tunneling between both leads. 
The clearly visible interference effects appear only for /il not too distant from Ed- 

In Fig. 6 we show the derivatives of the average current vs. the QD energy level Ed 
obtained for some values ol jiL- These are the results of the intersection of the surfaces 



given in Fig. 2 with the planes at constant values of /x^ or the results of the differentiation 
of curves shown in Fig. 4. Again, the most visible differences between the results obtained 
for Vrl — and Vrl ^ are present for the QD energy level Sd localized approximately 
between chemical potentials of both leads (compare, for example, the curves calculated 
for = 8). 

Fig. 7 presents the comparison of the d{jL{t)) / d^L vs. curves calculated for 
vanishing Vrl (left panels) and for Vrl — 10 (right panels) for two different values of 
the amplitudes (A^ = Ai/2, A^? = 0). At the vanishing value of Vrl, the shape of 
the curves is symmetrical in relation to the values fJ^L = although for greater A^ some 
shoulders appear on both sides of the corresponding peaks in the distance ^ Al/2 from 
the curve centres. For non- vanishing Vrl, the corresponding curves are approximately 
asymmetric and for large /il they tend to constant, non- zero values corresponding to 
linear increasing of the current at large /il- It is interesting that with the increasing 
amplitudes A^ and A^ very clear structures appear on both sides of the corresponding 
curves. Note, that all these curves can be obtained, for example, from the one calculated 
for Ed — O and moved along the /Xi-axis by the corresponding value (equal to Sd)- 

Fig. 8 shows the dependence of the average current {jL{t)) on the QD energy level 
and the direct coupling Vrl between both leads. The upper and bottom panels correspond 
to different values of the amplitude Al (A^ = Al/2). A very distinct transition from 
the symmetric to nearly antisymmetric behaviour of the current vs. ea curves is observed 
with the increasing value of the over-dot additional coupling between the leads for both 
values oi Al- The larger amplitudes of the left lead and QD levels oscillations result only 
in some broadening of the characteristic features of the average current vs. Vrl and £d 
surface and do not introduce any additional structures on these surfaces (for the range of 
parameters where cu ^ T and Hl — IJ'R^s not very small). 

Fig. 9 presents the average current (jl(^)) obtained for the case in which only the 
QD energy level Sd oscillates with some frequency (a; > F ) and for the small source- 
drain voltage Hl — IJ'R — 0.2 . The subsequent panels a, b, c and d correspond to the 
increasing value of the amplitude A^ and the different curves (broken, thin, thick and very 
thick) describe the QD without the over dot channel, Vrl = 0, and with this channel at 
Vrl = 4, 7 and 10, respectively. For vanishing Vrl (broken curves) we observe for small 
amphtude A^ only a central resonant peak (Fig. 9a). With increasing A^^, the subsequent 
peaks appear and the distance between them and the central peak is an integer multiple 



of the frequency cu (sidebands). The location of peaks is independent of the amphtude 
All but their relative intensity values change and with increasing A^; the height of the 
central peak is reduced. For the relatively large amplitude A^, the heights of the two 
neighbouring peaks (at Sd ± u) arc approximately equal to the height of the central peak. 
If we take the additional over-dot tunneling channel into consideration, especially for 
small Arf, the asymmetric shape of the current curve is observed and this asymmetry 
increases with the increasing strength of the over-dot coupling between both leads. With 
the increasing amplitude A^ this asymmetry is reduced largely due to the extra, photon- 
assisted tunneling peaks whose strength increases with the increasing A^. For sufficiently 
large values of A^ and Vrl the functional dependence of the average current on the QD 
energy level Sd is nearly the same for vanishing and non-vanishing over-dot tunneling 
channels. There is only one difference - for large Vrl the corresponding curve is shifted 
to the higher values due to the direct channel between both leads. 

The last two Figures 10, 11 are devoted to the analysis of the average current 
{jiit)) dependence on the oscillation period 27f/uj of the external fields. In Fig. 10 
we show the overall dependence of {jhit)) on 27: /uj and the QD energy level Sd- The 
upper (lower) part of the Figure presents the results for the vanishing (non-vanishing) 
over dot channel between both leads. The most visible differences between averaged 
currents calculated for Vrl = and Vrl 7^ can be observed for Sd < especially for 
small values oi2Ti/uj. More detailed analysis of the {jhit)) dependence on the oscillation 
period 2t: /uj is shown for some chosen values of the parameters Sd and Vrl in Fig. 11. 
The thin (thick) curves correspond to Sd — ^ {sd — 1) and the solid (broken) curves 
correspond to Vrl — 4 (Vrl — 0). Additionally, we give the resutls for two values of 
the amplitude A^ (A^^ = Al/2, Ar = 0), i.e. for A^^ = 5 and A^ = 10, the left and 
the right parts of Fig. 11 respectively. We observe the characteristic average current 
oscillations damped with increasing 27r/a; for ed lying in the central part between the left 
and right chemical potentials (see also [9] ) . These oscillations are present for both Vrl — 
and Vrl ^ and more visible for greater amplitudes A^ and A^ but the maxima and 
minima of the oscillating average current are localized at the same value of the oscillation 
period. Note, that the existence of the additional over-dot tunneling channel results 
approximately in shifting the corresponding curves to higher values without any additional 
serious modifications. For the QD energy level lying away from the middle point between 
the left and right chemical potentials, the average current is still an oscillating function of 



27r/a; although these oscillations are less transparent and their oscillation period is much 
longer. 

To conclude, we have provided a detailed investigation of a QD connected to two 
leads with an additional over-dot tunneling channel. The harmonic external microwave 
fields were considered as applied to the QD and two leads which result in time-dependence 
of the corresponding QD and leads energy levels. The QD charge and the average cur- 
rent flowing in this system were calculated within the evolution operator technique. The 
corresponding matrix elements of the evolution operator required for the calculation of 
the QD charge and current were presented in the form of the infinite series of multiple 
time integrals of the functions containing the information about the coupling between 
the QD and leads or in the form of the integro-differential equation. Applying the WBL 
aproximation we were able to obtain all required evolution operator matrix elements in 
closed forms and give the final analytical expressions for the QD charge and current. We 
have performed the extended numerical calculations for the QD charge, the average cur- 
rent and the derivatives of the current with respect to the gate and source-drain voltages 
(Fig. 4). The most spectacular infiuence of the additional bridge tunneling channel is 
visible in the {jL{t)) dependence vs. the position of the QD energy level at the constant 
source-drain voltage. Going from the vanishing values of the over-dot tunneling channel 
strength to the non-vanishing one due to the interference effects, the corresponding curve 
transforms from the symmetric to nearly antisymmetric shape. Similar influence of the 
non- vanishing Vrl is visible in the dependence of d{jL{t)) /d/iL vs. /il (Fig- 7). The 
characteristic behaviour of the average current vs. the QD energy level position at the 
small source-drain voltage is observed for the case when the external oscillating field is 
applied only to the QD (Fig. 9). For the vanishing over-dot tunneling channel at the small 
amplitude A,^, the main resonant peak is only observed and with the increasing amplitude 
Ad the next peaks localized at Ed equal to the multiplicity of cu appear corresponding to 
the photon- assisted tunneling. In that case the current is a symmetric curve centered 
around the main resonant peak. For the nonvanishing over-dot tunneling channel and 
small amplitudes Ad, the character of the dependence of the average current on the QD 
energy level position transforms with increasing Vrl from the symmetrical to nearly anti- 
symmetrical behaviour. This tendency to the antisymmetrical behaviour with increasing 
Vrl at small amplitude A^; is reduced with increasing Ad- For sufficiently large amplitude 
Ad the overall behaviour of the average current vs. Ed is very similar for different values 



of Vrl and does not manifest the tendency for the antisymmetry with increasing Vrl- 
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Figure 1: The average current {jL{t)) against the left lead chemical potential fj,L and QD 
energy level ea for Vrl — (the QD without the over-dot channel, the upper panel) and 
for Vrl = 10 (lower panel), /ir = 0, V = 4, = 2, Aa = 1, Ar = 0, uj ^ 2 and all 
energies are given in F units. 





Figure 2: The derivatives of the average current with respect to the QD energy level e^i 
d{jL{t)) /dsd, against /il and Sd for Vrl — (upper panel) and for Vrl — 10 (lower panel). 
The other parameters as in Fig. 1. 





Figure 3: The derivatives of the average current d{jL{t)) /d/iL, against fj,L and Sd for 
Vrl = (upper panel) and for Vrl = 10 (lower panel). The other parameters as in Fig. 1. 
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Figure 4: The average current {jhit)) against for given values of (beginning from 
jiL — —4 up to Hl — The left and right panels correspond to Vrl — and Vrl — 10, 
respectively, and the other parameters as in Fig. 1. 




Figure 5: The average current (jL(i)) against /il for given values Ed (beginning from 
Ed — —4 up to Ed — S). The broken curves correspond to Ed — 0. The left (right) panel 
corresponds to Vrl — {Vrl — 10). The other parameters as in Fig. 1. 
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Figure 6: The derivatives of the average current d{jL(t)) /dsd with respect to the QD 
energy level for given values of fi^ (beginning from fi^ = —4 up to /i^ = 8). The 
broken curves correspond to hl = —4. The left (right) panel corresponds to Vrl = 
{Vrl — 10) and the other parameters as in Fig. 1. 








Figure 7: The derivatives of the average current with respect to fi^, d{jL{t))/ dfi^, for given 
values of Ed (beginning from Ed = —4 up to = 8). The left (right) panels correspond 
to Vrl — {Vrl — 10) and upper (lower) panels correspond to A^, = 2 (A^, = 4). The 
other parameters as in Fig. 1. 
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Figure 8: The average current {jL{t)) against Vrl and ea- The upper (lower) panel 
corresponds to = 2, = 1 (Al = 4, A^ = 2). = 2, V ^ 4. 




Figure 9: The average current {jL{t)) against Ed for the oscillating QD energy level at 
Vrl — 0, 4, 7, 10 - broken, thin, thick and very thick curves, respectively The panels a, 

b, c and d correspond to = 1, 3, 5 and 7, respectively, uj = 5, T = 1, V = 4, hl = 0.2, 
Al^Ar^ 0. 
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Figure 10: The average current {ibit)) against Sa and the period of the time oscillation 
of the external field for Vrl = - the upper panel and for Vrl = 10 - the lower panel, 
y = 4, = 10, = 10, Arf = 5, Afl = 0. 



Figure 11: The average current {ibit)) against the period of the time oscillation of the 
external field for Vrl = - the broken curves, and for Vrl = 4 - the solid curves. The 
thin (thick) curves correspond to 8^ = ^ {sd = 5). The left (right) panel corresponds to 
Al = 5, Ad = 2.5 {Al = 10, Ad = 5). Ar^ 0, V ^ ^ 10. 



